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Abstract — Brain MRI images are mostly used in order to 
identify the diseases related to brain. Of late many models are 
presented by utilizing the MRI medical imaging technology. 
However most of the models failed to identify the inhomoginities 
in the brain accurately. Among the various inhomoginities, 
seizures are responsible for alzimers diseases which cause brain 
deaths. Hence this paper addresses a mechanism for the early 
identification of seizures by reducing the dimensionality of the 
images using PC A. The database from UCI is considered for 
experimentation and the performance is evaluated using 
benchmark metrics like image fidelity, mean square error, peak 
signal noise ratio (PSNR). 


Index Terms — Brain images, Alzimers disease, seizures, 
performance evaluation, PCA, Quality Metrics. 


I. Introduction 

Brain imaging is a mainly focused area at the present hour 
due to the increase in mortality rate due to the brain related 
diseases. Many models are proposed by the researches to 
identify the diseases related to brain. Among the various 
diseases in brain, multiple Seizures, lesions, Selorisis and 
inhomoginities corrections are mainly focused areas. The 
brain image consists of three vital tissues namely White 
Matter (WM), Grey Matter (GM) and cerebrospinal fluid 
(CSF). The MS Seizures is mainly witnessed in young adults 
and it affects the white matter. MS lesions affect the grey 
matter part and seizures are the outer layers in any brain 
images. In order to have a complete study of the brain and 
identify the related diseases, it is customary to extract the 
images accurately and the inhomoginities in the intensity of 
the pixels lead to the mis-interpretations. Hence, affective 
mechanisms are to be developed for the identification of these 
diseases. This article in particular highlights the contribution 
made in the identification of seizures. In any brain MRI 
images, tl and t2 weighted images are mostly witnessed. The 
high intensity levels inside these images reflect the white 
matter (WM), the medium intensity levels reflect the grey 
matter (GM) and the lesser intensity levels help towards the 
identification of cerebrospinal fluids (CSF). Any intensity 
levels which differs the above levels are formulated as the 
outer layers or the abnormal pixel intensity inside the medical 
image regions are coined as seizures. Many models have been 
proposed by the researchers to witness the seizures. An 
identification failure of the seizures accurately leads to 

Manuscript received January 19, 2015. 

M Vamsi Krishna, Research Scholar, CUTM. 

Dr. P. V. G. D. Prasad Reddy, Professor, Department of CSSE, Andhra 
University, prasadreddy 

Dr. Y. Srinivas, Professor, Department of IT, Gitam University, sriteja 

Dr.S.S.Nayak, Dean R & D, Centurian University 


abnormalities in brain or leads towards specific diseases like 
alzimers. Hence it is customary to identify the diseases more 
accurately. The global threshold of the medical brain image is 
extracted and all the pixels above this global threshold are 
generally categorised as seizures. Good amount of research is 
projected to identify the seizures based on intensity levels 
using gaussian mixture models (GMM). However, if the 
image is not pre-processed and contains noisy data, it is 
difficult to sensitize the noisy part and the abnormal intensity 
value which is supposed to be the seizure pixel. Hence the 
shape and density of the pixel should be considered along 
with the intensity. Therefore this article utilizes these features 
and proposes a model based on bivariate gaussian mixture 
model. The dimensionality reduction technology of Principal 
Component Analysis (PCA) is considered for reducing the 
dimensionality and extraction of the feature vectors. The rest 
of the paper is organized as follows. Sections 2 of the paper 
deals with bivariate gaussian mixture model (GMM), the 
features considered are presented in section 3. The 
dimensionality reduction using PCA is proposed in section 4. 
Section 5 of the paper highlights the experimentation together 
with the results obtained. Section 6 of the paper highlights the 
performance of the proposed model and section 7 of the paper 
summarizes the paper. 

II. BIVARIATE GAUSSIAN MIXTURE MODEL 

In order to have an effective recognition rate, in identifying a 
tumour, one need to consider a bivariate approach since it can 
handle the characteristics of the tumours more effectively. 
Any tumour processes a bivariate approach rather than a 
univariate approach. 

The probability density function of the bivariate gaussian 
mixture model is given by 

f{xl,x2,p) = [27t v '{I-p2) - 1 e (^^)(xU-2{ixlx2 + x22)] 

Where xl , x2 are the parameters of the tumours, {.t is the mean 
of the MRI image under consideration, p is the shape 
parameter. For a particular value ofp — 0, it assumes the 
circular normal distribution and for p = 0, and ffl = tr2, 
it formulates an elliptical normal distribution. 

III. FEATURE EXTRACTION 

In order to demonstrate the proposed methodology, a 
database is extracted from UCI brain imaging which contains 
the information pertaining to various brain images with 
diseases. The image pixels are categorized into four groups, 
where the outer group consists of the pixels having very high 
intensity values or abnormal pixels. The next intensity levels 
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are assumed to be the white matter (WM), grey matter (GM) 
and cerebrospinal fluid (CSF) respectively. Each pixel in the 
brain image under consideration are extracted and given as 
input to the bivariate gaussian mixture model (GMM) 
discussed in the section 2 of the paper. The features are the 
pixel intensity levels and the Eigen values obtained after using 
PCA. 

IV. PRINCIPAL COMPONENT ANALYSIS (PCA) 

PCA is utilized in this paper to reduce the dimensionality of 
the image and identify the principal components called the 
Eigen values. Here each image is considered are categorised 
into four regions namely region of cerebrospinal fluid (CSF), 
grey matter (GM), white matter (WM) and outer layer. This 
region separation is highlighted by plotting the histogram of 
the image or from the frequency cure of the image. Each 
image consists of a mixture of sub images each exhibiting a 
particular distribution. Hence the whole image is assumed to 
be a mixture of distribution. The PCA is applied to each of 
these regions and the Eigen values along with the intensity 
values are given as input to the bivariate model given in the 
section 2. 


V. METHODOLOGY 

In order to evaluate the proposed method, the images from the 
UCI dataset are extracted and are normalised to a size of 100 
by 100. Each image is pre-processed for identifying the 
missing data and eliminating the noise. These processed 
images are extracted to identify the regions. On each of these 
regions, along with the pixel intensities of the region, the 
probability density function of the bivariate gaussian mixture 
model is generated. Any pixel is categorised basing on the 
likelihood estimates. If the maximum likelihood of a pixel 
satisfies a region, the pixel belongs to that region. This 
mechanism highlights the seizures more accurately because 
we are considering the likelihood of each pixel against each 
region more abruptly. The intensities of the voxel pixels are 
the outer layers can be easily identified by this method. 

A. SEGMENTATION ALGORITHM 

After estimating the initial parameters, the first step is 
assigning the pixels to the segments. The procedure is carried 
out by the segmentation algorithm, consisting of the following 
steps. 

Step 1: Acquire the pixel intensities of the image. Denote 
them by Xy. 

Step 2: Find initial basing on the intensity values. 

Step 3 : For each region obtain the initial estimates for mean 
Pi and variance Gi. 

Step 4: Implement the segmentation algorithm by 

reconstruct the image by assigning the images based on 
maximum Likelihood estimate. 

Step 5 : Estimate the image quality by using metrics from the 
reconstructed image using the step 4. 


VI. EXPERIMENTAL RESULTS AND PERFORMANCE 
EVALUATION 

After developing the segmentation algorithm, the algorithm is 
applied to T1 -weighted brain medical images obtained from 


the UCI database of dimensions 150x174 and 163x199 
respectively. The initial parameters obtained are used for the 
reconstruction by assigning each pixel in the PDF .In this 
paper a methodology for Seizures identification is conducted 
basing on BGMM. The brain image is classified into three 
major regions white matter (WM), grey matter (GM) and 
cerebrospinal fluid (CSF) and the likelihood estimates of each 
of these regions are estimated. Each pixel is compared against 
the region and if it does not satisfy the likelihood estimates, 
they are classified as a Seizures. This process is repeated for 
all the images. In order to evaluate the performance of the 
reconstructed image, image quality metrics are used and the 
metrics utilized for this purpose are presented below. 


Table 1. Quality Measures 
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From the above Table 2, it can be clearly seen that the model 
developed by using K-Means clustering shows better results 
with respect to the quality metrics. The model is compared the 
existing models based on Gaussian Mixture Model 
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VII. PERFORMANCE EVALUATION 

In order to evaluate the above model, quality metrics 
like image fidelity, peak signal noise ratio, mean squared 
error, average distance are used. The formulae for calculation 
of these metrics are given as follows 


Quality 

Metrics 

Formulae to Evaluate 

Average 

Difference 

M jY 
j= 1 iFl 

Where M,N are image matnx row's and columns 

Maximum 

Distance 

rrm[|F(/4) -F(j,ii)|] 

Image 

Fidelity 
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VIII. SUMMARY 

A medical image segmentation technique based on Bivariate 
Gaussian distribution model with K-Means clustering is 
developed and evaluated. The results obtained by this 
algorithm outperform the existing method of GMM. This 
method can be mainly suited in particular cases of medical 
pathology for the identification of SEIZURES accurately 
there by helping in proper diagnosis and preventing 
disabilities. 
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